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QQ ■ We investigate the long time behaviour of the one-dimensional ballistic aggregation model that 

' represents a sticky gas of A'^ particles with random initial positions and velocities, moving determin- 

, istically, and forming aggregates when they collide. We obtain a closed formula for the stationary 

^SJ ■ measure of the system which allows us to analyze some remarkable features of the final 'fan' state. 

In particular, we identify universal properties which are independent of the initial position and 
velocity distributions of the particles. We study cluster distributions and derive exact results for 
) . extreme value statistics (because of correlations these distributions do not belong to the Gumbel- 

• Frechet-WeibuU universality classes). We also derive the energy distribution in the final state. This 

' model generates dynamically many different scales and can be viewed as one of the simplest exactly 

solvable model of A'^-body dissipative dynamics. 

^ ■ PACS numbers: 68.43.Jk, 02.50.-r, 05.40.-a, 47.70.Nd 
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I. INTRODUCTION 



^ , The ballistic aggregation model is one of the simplest interacting particles processes of non-equilibrium statistical 

■ mechanics and as such has attracted a lot of attention in the past decades This model represents a gas of N 
unit-mass particles forming clusters through perfectly inelastic adhesive collisions. The motion of a particle between 
two collisions is deterministic and free (i.e. ballistic) and the total mass and momentum in a collision are conserved 
whereas the kinetic energy is dissipated. The stochasticity in this model is due only to the initial configuration 

, which consists of single particles randomly located with uncorrelated random velocities drawn from a continuous 
(~| ■ distribution. This dissipative system, usually referred to as ballistic aggregation or sticky gas, appears as a minimal 
O '. model of cluster formation and provides a relevant statistical description of the merger of coherent structures such as 
vortices, thermal plumes, flowing granular media [3] or the accumulation of cosmic dust into planetoids. The ballistic 

■ aggregation model also plays a role in the study of the large-scale structure of the universe because the motion of 
fS| self-gravitating matter in the expanding universe is similar to that of matter moving solely by inertia 0, jjj] . Another 

^ ' noteworthy feature of this model stems from its connection with the Burgers equation, an important toy model in 
OO , the study of turbulence: it can be shown that at very high Reynolds number, the solution of the Burgers equation 

■ in the long time limit consists of a series of shocks which follow exactly the dynamics of the ballistic aggregation 

g, model gag. 

Ballistic aggregation exhibits at long time t a self-similar coarsening behaviour first studied by Carnevale, Pomeau 
and Yoimg In one dimension, the coarsening regime occurs for intermediate times satisfying t ~ TV^^^, N being 
the total mass or equivalently the total number of initial particles in the system. For longer times, and in the absence 
of any boundary, the system ultimately reaches a state where no more collisions are possible: the particles are grouped 
into clusters of different sizes (or masses) with velocities increasing from left to right. This final ordered state will be 
called the "fan" state. Once this state is reached, there is no further loss of energy and the system becomes stationary. 
The aim of this paper is to investigate the properties of this fan state. We find an exact analytic formula for 
; ^ ' the joint probability distribution of the sizes and velocities of the clusters by mapping the final state of the ballistic 
aggregation model to the convex minorant of a one-dimensional random walk. From this invariant measure that fully 
characterizes the set of all possible fan states, various statistical properties of the model in the long time limit will be 
derived. In particular, we shall retrieve the known fact [ill [l^ that the probability of obtaining a fan state with 
exactly k clusters is a purely combinatorial factor (which is universal in the sense that it does not depend upon the 
initial velocity distribution of the particles). We shall prove that the cluster distribution in the fan state is identical 
to the statistics of cycles in a random permutation of N objects [isl. [l^L both problems having a common underlying 
combinatorial structure related to the convex minorant construction This will allow us to characterize the 

different scales that are dynamically generated in the system: in the large N limit, the size of a typical cluster is of 



order N/\nN, the largest cluster contains a finite fraction of the total mass N and hence grows linearly with N, the 
smallest cluster scales as IniV, and the rightmost cluster, which has the leading velocity, has a size of order ^/N. We 
shall also derive the (non-universal) joint distribution of mass and velocity of this 'leader' which is very different from 
the usual extreme-value distributions obtained for uncorrelated random variables. Finally, we shall show that the 
energy remaining in the system after all collisions have ended is of the order In N and shall calculate the distribution 
of the scaling variable E/\nN. 

The paper is organized as follows. In Section-II we define the ballistic aggregation model and focus on its final 
fan state. In particular, we show how the statistical properties in the fan state are related to those of the convex 
minorant of an associated one dimensional random walk problem. This mapping allows us to derive explicitly the 
joint distribution of the cluster sizes and their velocities in the fan state which is one of the key results of this paper. 
In Section-Ill we show how this joint distribution can be used to derive the cluster size distribution in the fan state 
which turns out to be universal (with respect to the initial velocity distribution). This universal property is shown 
to be related to the statistics of cycles lengths in a random permutation problem. In Section-IV, we further exploit 
the joint distribution to compute exact asymptotic results for certain extreme variables in the fan state such as the 
size and the velocity distribution of the rightmost ('leader') cluster. In Section-V, we derive the exact asymptotic 
distribution of energy in the fan state. Finally we conclude in Section- VI with a summary and some open problems. 

II. STATISTICAL DESCRIPTION OF THE STATIONARY STATE 
A. Definition of tiie Model 

The ballistic aggregation problem in one dimension is an elementary interactive particles process, the dynamics of 
which is defined as follows. At the initial time t ~ 0, we start with N particles randomly located on an infinite line. 
We label the particles as i = 1,2,. ..,iV according to the order of their initial positions from left to right (see Fig. 1 
where we have taken = 6). We denote by Vi and rrii, respectively, the initial velocity and the mass of the i-th 
particle. For t > 0, each particle moves ballistically at constant speed. When two particles meet they coalesce and 
form a single new particle whose mass is the sum of the initial masses and whose momentum is the sum of the initial 
momenta. The only randomness in the model lies in the initial conditions, i.e. the initial positions and velocities 
of each particle. We shall assume that the initial velocities {wi}'s are uncorrelated random variables, drawn from a 
common continuous probability density function (PDF) (f>{v). For simplicity, we assume that all the initial masses 
are the same and take them to be unity. Once the initial conditions are given, the evolution of a given realization is 
deterministic and fully determined by the mass and momentum conservation laws. The shocks being totally inelastic, 
energy is dissipated at each collision. 

At late times t 3> 0(1) (but before the system feels the finiteness of the mass N), the ballistic aggregation model 
exhibits a coarsening regime where typical cluster size grows with time as a power law, as studied first by Carnevale, 
Pomeau and Young 1] using scaling arguments and numerical simulations. For example, in this regime, the average 
size (or mass) of a particle-cluster grows as ^2-0/(2-1-13) velocity decreases as jj being the spatial 

dimension. The global mass of the system being conserved, the total energy of the system decays as (2+^') _ xhese 
exponents derived in [l| turn out to be correct in one dimension where the problem is exactly solvable [1,0, This, 
however, turns out to be a bit fortuitous. In higher dimensions, the exponents predicted in_Jl| turn out to be incorrect 
due to strong correlations between the velocities of the colliding clusters at late time s |l7| . In one dimension, it is 
possible to calculate exactly (in one dimension) the mass distribution of the clusters 0, H 0] in the scaling limit when 
TV — > oo {N being the total mass of the system) and i — > oo, keeping the ratio t/N^''^ finite 0, 0, However, these 
scaling results for the coarsening regime are valid only for intermediate times satisfying t ~ N^^^. When t ^ N^^^, 
the system evolves into a stationary state in which no more collisions can occur: the particles are grouped in k disjoint 
clusters of different masses, where k itself is a stochastic variable. Each cluster moves at a constant velocity, and 
the speed of a given cluster is larger than that of its left neighbour (if any) and less than that of its right neighbour 
(if any). In this ultimate state the clusters keep on moving farther apart, i.e., they fan out from each other with 
increasing time, thus justifying the name "fan" state (see Fig. 1). Once this fan state is reached, there is no further 
collision and hence no further dissipation of energy. 

In this work, we shall focus on the statistical description of the fan state of the random aggregation process. We 
shall show that properties of the fan state related to the sizes of the clusters (regardless of their respective ordering and 
their velocities) are universal as they do not depend on the PDF (j){v) from which the initial velocities of the particles 
are drawn. The cluster statistics can in fact be mapped to the cycle length distribution in random permutations, 
which are fundamental combinatorial objects. From this observation the distribution of the cluster masses and in 
particular the typical mass of the largest and the smallest cluster can readily be calculated. When the distribution 
of the velocities of the clusters is taken into account, universality with respect to 4>{v) is lost. However, in the large 
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FIG. 1: An example of a configuration of 6 particles with identical initial mass rrii = 1 and varying initial velocities, evolving 
in the long time limit into the fan state: the velocities are now increasing from left to right and there are no more coUhsions. 

N limit, we will see that some universality is restored, thanks to the central hmit theorem, in the size distribution of 
the rightmost cluster (the leader) provided the second moment cr^ = (f>{v) dv is finite. Finally, we are also able 

to calculate explicitly the energy distribution in the fan state, but only for the special case when (t>(v) is Gaussian. 

B. The Fan State as a Convex Minorant 

The fan state of the random aggregation process can be determined geometrically from the initial conditions 
We now explain this graphical interpretation of the fan state as a convex minorant of a random walk. Recall that we 
label the particles as? = l,2,...,iV according to the order of their initial positions from left to right (see Fig. 1) and 
that we denote by Vi respectively the initial velocity of the i-th particle (all the initial masses being unity). For any 
realization of the initial state {f i, V2, ■ ■ ■ , vn}, the final state is uniquely determined i.e. the number of final clusters, 
their masses and their velocities are unique. 

The initial state {vi,V2, ■ ■ ■ , vn} of the system is represented by a broken graph {Pq, Pi, ... , Pn) (see Fig. 2) such 
that: 

(i) Po = (0, 0) is the origin; 

(ii) the coordinates of Pi are given recursively by Pi = Pi-i + (1, Vi) for 1 < i < TV. 

In other words, the speed Vi of the i-th particle is represented by the slope of the line (Pi^iPi). We also emphasize 
that the horizontal coordinate of Pi does not correspond to the actual position of the i-th particle but only on its 
label. 

More generally, if we had started with particles having different masses {mi,TO2, • • ■ i'tiat} the coordinates of Pi 
would be Pi = Pi-i + mi{l,Vi) [Tl|. In other words, the initial state is drawn as a random walk in the cumulative 
momentum and cumulative mass space as shown in Figs. [5] and [31 

Suppose that the first collision occurs between particles i and {i + 1) with velocities Vi and Vi^i respectively. 
These two particles can aggregate if Vi > Wi+i, which means that the slope of the [Pi-iPi) is larger than the slope 
of the {PiPi+i). Equivalently, this means that Pi is located above the segment (Pi_iPi_|_i) i.e., locally, the graph 
{Pi^iPiPi+i) has a negative curvature. After this collision, a cluster is formed with mass 2 and velocity {vi -\-Vi+i)/2. 
This cluster is represented by the vector {Pi-iPi+i) (see Fig. [2]) with coordinates (2,Di + Wi+i). We note that the 
slope of this vector again represents the velocity of the cluster. We now have a system of TV — 1 'particles', with N — 2 
particles of mass 1 and one particle of mass 2. The state of this system is represented by a broken line in which the 
angle (Pj_iPiPi+i) with downward curvature is replaced by its base (see Fig. [2]). Similarly, the next collision is also 
represented graphically by replacing another angle with downward curvature by its base (see Fig. [2]) . This process 
will continue iteratively and the particles will aggregate forming clusters till all angles with downward curvature have 
been eliminated, i.e. all collisions have occurred. It follows that for any given initial state, the final state is uniquely 
given by the convex minorant of the corresponding random walk. Each line segment of the convex minorant represents 
a cluster in the fan state, — the horizontal and the vertical components of the segment give respectively the mass and 
the momentum of the cluster, i.e., the slope gives the velocity of the cluster. The geometry of the convex minorant 
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FIG. 2; Initial configuration is represented by a graph (thin line) joining Pq to Pjv (we have taken rrij — 1 Vj). This graph can 
be interpreted as a random walk in the cumulative momentum and cumulative mass space. After the first collision the particles 
i and i + 1 aggregate into one cluster: in the associated graph, the angle (Pi-iPiPi+i) is replaced by the segment (Pi-iPi+i). 
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FIG. 3: We draw the convex minorant of the graph representing the same initial configuration as in Fig. [5] This convex 
minorant represents the final 'fan' state in which no more collisions can occur. Each straight segment of the convex minorant 
corresponds to a cluster of n; particles with momentum niUi. 

(including the number of segments which represents the number of final clusters), however, differs from one realization 
of initial conditions to another. Various statistical properties of the ballistic aggregation in the fan state can be found 
from the statistical properties of the underlying convex minorant. 

C. The Invariant Measure in the Fan State 

We now derive the joint probability distribution which describes the fan state. Let pAr(fc; {ui, ui}) denote the joint 
PDF of having k clusters in the fan state of masses and velocities rii and Ui respectively, with i = 1, 2, . . . , fc. Note that 
the various segments of the convex minorant (see Fig. are correlated as their slopes must be increasingly ordered 
from left to right, i.e., ui < U2 < • • • < Uk- Moreover, the total mass is conserved, X^iLi ~ ^ ■ Once these two 
constraints are specified, the momentum of each final cluster is essentially determined by the conservation law: 

riiue = ^ rrijVj (1) 

for the ^-th cluster with £ = 1,2, ... ,k. We have defined here the cumulated masses Nk = ni + . . . + Uk- The sum in 
the above equation is thus restricted to the initial velocities of those particles which form the l-th cluster. We also 
recall that rrij = 1 Vj. However, most importantly, to form a specific (say ^-th) final cluster (i.e., a line segment in the 




underlying convex minorant) by aggregation of ng consecutive initial particles, the portion of the underlying random 
walk of Hi steps must be such that it originates from one end of the line segment of convex minorant and finishes at 
its other end, while staying above the segment in the intermediate steps. Now, consider a random walk of ni steps 
that starts from one end of a line segment and finishes at its other end, but is otherwise allowed to cross the segment 
in the intermediate steps. For any realization of such a walk, given that there is one unique minimum with respect 
to the segment, if we consider all the ng cyclic permutations of the steps, then out of ng different realizations of the 
walks there is one and only one arrangement where in the intermediate steps the walk stays above the segment. This 
property is known as Raney's lemma [isl . [l9j . (Here, the uniqueness of the arrangement is guaranteed by the fact 
that the PDF (j){v) is continuous). 

Thus, the probability that a random walk of ng steps that starts at one end of a line segment and finishes at its 
other end while staying above the segment in-between equals l/ng. Note that this probability is independent of the 
slope ug of the segment as well as the PDF (p{v) of the jump-lengths. Finally, gathering all the above inputs together 
we write 

Cfc \ k-l 
i=l I 1=1 

where the angle brackets (• • • ) denote the averaging over the set of initial velocities {wi}. This key result provides a 
full statistical description of the clusters sizes and velocities in the fan state. From this joint probability distribution, 
all the statistical properties of the final state can be derived. 



III. UNIVERSAL STATISTICAL PROPERTIES OF CLUSTERS IN THE FAN STATE 

In this section, we compute the distribution of the number and sizes of the clusters in the fan state: we will show 
that this distribution, obtained by integrating out the velocities of the clusters from the joint probability distribution 
given in Eq. is universal with respect to the initial continuous PDF (^(v). 



A. Statistics of the Total Number of Clusters in the Fan State 



We first calculate pjv(fc), the probability that the final state contains k distinct clusters. The value of pjv(fc) is 
obtained by summing Eq. over all ni > O's and integrating it over m^'s, keeping fc fixed. Thus, 



fc-i fc 
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In this equation, we have used 
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where (j) is the distribution of the initial velocities. We note that the function V{n, u) is normalized: 
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duV{n,u)^l. (5) 
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To disentangle the constraints in Eq. ([3]), it is useful to consider the following generating function: 

N {m} i=l i=l i=l i=l 



with 



Pziui) = V{ni,Ui) . 

^ — ' Uj 



The integral on the right hand side (r.h.s) of Eq. ^ can be evaluated recursively and is found to be 
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where the last equality is obtained by using the definition of pz (u) in ([7]) and the normalization property in 
for any given value of fc > 0, we have 

1 
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After extracting the coefficient of z on the r.h.s of this equation, we conclude that 
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(8) 
Thus, 
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This formula can be expressed in terms of some classical combinatorial numbers as follows. It is convenient to 
introduce a second auxiliary variable t and to calculate the generating function of Eq. ([9|) with respect to t: 



k N k ' 



t'^ (1 - z)-^ = l + tz + 



t{t + l) 2 t{t + l){t + 2) 3 



2! 



-z^ + 



3! 



Z-' + 



Identifying the coefficients of z on both sides of this equation we obtain the identity: 
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Using the classical formula 
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t{t + l)(t + 2)...{t + N -l)=^Si{N,k)t^ 



k=l 



(13) 



that defines the unsigned Stirling numbers of the first kind 20] , Si {N, k) , we conclude that the probability of finding 
k clusters in the fan state for a system starting with N distinct unit-mass particles is given by 



PNik) = = 



Nl 



with fc = 1,2, ... ,7V. 



(14) 



We thus recover by a different method the result of Ref . [ll| . 

From Eq. (fT3|l . we observe that X^fcLi '5'i(^: fc) = Thus, p]\f{k), and consequently the joint PDF given by 
Eq. are both clearly normalized as it must. Furthermore, using Eq. p^ . one deduces that the mean number of 
clusters is given by 1 + ^ + . . . + . 

In fact, it is well known 20] that the unsigned Stirling numbers Si{N, k) enumerate the number of permutations of 
TV elements with k disjoint cycles exactly: thus, pN{k) is identical to the probability distribution of number of cycles 
in random permutations with uniform measure. Many results have been derived concerning the statistics of random 
permutations. For example, it has been shown in the context of permutation cycles [13] that for large N, pN^k) has 
a Gaussian form around its mean (fc) ~ In A^, with a variance which is also of order In A*". Therefore, the typical mass 
of a cluster is of order N/\nN. Using this together with the t^^^ asymptotic growth of cluster sizes before the fan 
state is reached, the characteristic time to reach this state is estimated to be tc ~ {N/ hiN)'^^^. 



B. Distribution of Clusters Sizes including the Biggest and the Smallest 



We now consider a more refined observable. It may be tempting to calculate the cluster size distribution i.e. to 
integrate out the velocities from Eq. ^ and determine the marginal PDF pN{k; {ni})- Unfortunately, this marginal 
distribution is not universal: it depends on the choice of the initial PDF (this fact can be verified by working 
out explicit examples). However, if one considers the statistics of the clusters according to their sizes regardless of 
their spatial ordering then the result is again (/)(w)-independent. 

More precisely, let cj for j — 1,2,...,N denote the number of clusters of size j in the fan state. By mass 
conservation, we have J2j J'^j — ^- Then, it can be shown that the joint probability, Pr{ci, C2, . . . , cjv} is given by 



N ^1 

Pr{ci, C2, . . . , c^} = <5(7V - ^ jc,) n -— T • (15) 

The derivation of this formula from the invariant measure in the fan state, Eq. ([5]), is given in the Appendix A. This 
expression is identical to the distribution of cycles in random permutations where Cj's are to be identified with the 
number of cycles of length j [2^ . 

The generating function associated with this probability distribution is defined as follows: 

GN{h,...,tM)^ Mci,C2,...,CN}tl' ...t"^ . (16) 

ci,...,cjv 

This function Gn can itself be embeded in a 'grand-canonical' generating function given by 

OO J, 

r(z, {tk}) - ^ G^(ti, . . . , t^)^^ = exp ( ^ ^) , (17) 

N fc=l 

where the last equality is obtained by using Eq. (jlSp . This expression embodies all the necessary information and will 
be useful for further calculations. 

For example, we can study the extreme cluster sizes. Let Qn{L) be the probability that the largest cluster size 
(mass) ninax < L and Rn{S) the probability that the smallest cluster size rimin > S. It follows from these definitions 
that the probability that the largest cluster is of size L is given by 

Pr(n,,ax = i) = gAr(i)-Qw(i-l), (18) 

and similarly 

Pr {Unun = 5) = Rn{S) ^Qn{S+1). (19) 

Using Eq. (fT5|) we find 

Qn{L)^ J2 Pr{ci,C2,...,CL,0,...,0} = GAr(ti = l,...,iL = l,0,...,0), (20) 

Ci,...,Cl 

and similarly 

Rn{S)^ Pr{0,...,0,cs,...,cw} = Gw(0,...,0,is = l,...,ijv = 1). (21) 

cs,...,cjv 

Using the grand-canonical generating function defined in Eq. (jl7p . we readily obtain 

QL(^)-5]Qjv(i)^'^-exp(^^J and i?i(z) ^ i?jv(5)z^ = exp ( ^ J . (22) 

AT \n=l / N \n=S ^ ) 

The mean sizes of the largest and the smallest clusters in the fan state containing N particles are given by 

(nmax)jv = L\Qn{L) - Qn{L - 1)] = ^[1 - Qn{L)] (23) 

L L 

and (nniin)jv = YS[Rn{S)-Rn{S+1)]=YRn{S). (24) 



For large N, the asymptotic behaviour of (rimax) n and (rimin) n can be obtained by analyzing the associated generating 
functions in the z ^ 1 limit. For example, using Eqs. (|22p we obtain 



N L \n=l / L \ n=L + l / 



(25) 



When z 1, the sums on the right-hand sides can be transformed into integrals and the desired asymptotics are 
obtained by extracting the leading singularity at z = 1. Thus, for large N, we obtain 

(nmax) ~ with A = exp (^-x -J ^ dy^ dx = 0.6243299885 . . . (26) 

This constant is known as the Golomb-Dickman Constant [llj. Similarly, We find, (rimin) ^ e~^lnN, where C = 
0.57721566 ... is Euler's constant, and we again note that, (rimin) and (rimax) are identical to the mean lengths of 
the shortest and the longest cycles, respectively, of random permutations and were already calculated in this context 
[fl [2^ . Recently, the same constants also appeared in the context of the statistics of longest and shortest lasting 
records of independent and identically distributed random variables [H, HJ] 



IV. LEADER STATISTICS 



In this section we compute the statistics of the size and the velocity of the 'rightmost cluster' in the fan state. 
This rightmost cluster will be referred to as the 'leader' since it moves with the highest velocity. Let Lpf{n, v) denote 
the joint PDF of the leader's size n and velocity v in the fan state. To obtain this PDF, we start from the basic 
microscopic joint PDF of all cluster sizes and their velocities in Eq. ^ and integrate out the sizes and velocities of 
all clusters except the leader. As usual, it is convenient to consider the generating function which then reads 



N ^ k=l '' 
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(27) 



where Pz{u) and 'P{n,v) are defined respectively in Eqs. ([T]) and The (fc — l)-fold integral over the u^'s can be 
recursively evaluated as in Eq. ([8]) and we get 
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Substituting this in Eq. (I27p and summing over fc we get a somewhat compact expression 
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For the subsequent asymptotic analysis, it turns out to be convenient to rewrite the expression in Eq. ([29)1 in a slightly 
different form. Using the definition of pz {u) from Eq. (jT)) and the normalization in Eq. ([SJ , it follows that 



Pz{u)du^Y.T / ^(^'"^ 
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du 



ln(l - z) 



Pz{u)du. 



Using this in Eq. 



gives 



N 



n{l — z) 



V^n, v) exp 



Pz{u)du 



(30) 



(31) 



From the joint PDF LN{n,v) of the leader's size and velocity in Eq. ([3T|l one can then obtain the marginals: 
CAr(n) = L]y{n,v)dv for the size PDF and Dm{v) = ^"^^i L]\[{n,v) for the velocity PDF of the leader. It is 
evident from Eq. (|31|) that LN(n,v), for any finite N, depends explicitly on the initial velocity distribution (f>{v) 
since both V{n,v) as well as Pz(m) depends on (j){v). This is unlike the distribution of the cluster sizes as derived 
in Section III which is completely universal, i.e., independent of 4>{v) for any finite N. So, a natural question is: to 



what extent the leader size and velocity distributions are universal, i.e., independent of the details of 4>{v)l We will 
see below that the universal property holds for the leader's size distribution, but only in the large N limit and for 
continuous and symmetric (j^iv) with a finite variance cr^ = v'^(j){v)dv. In that case, for large N, one can use 
the central limit theorem and the limiting marginal size distribution C'N{n) essentially becomes universal, though the 
marginal velocity distribution Dn{v) still remains nonuniversal even in the large N limit. A natural question is how 
these results get modified when (f){v) is such that its variance is infinite. A particular example of this class is the 
Cauchy distribution for We will show that the Cauchy case is exactly solvable and one can derive the size and 

the velocity distribution of the leader explicitly. In subsections A and B below, we consider, for finite cr^, the limiting 
size and the velocity distribution of the leader respectively. In subsection C, we derive the explicit distributions for 
the Cauchy case. 



The Limiting Leader Size Distribution 



In this subsection we compute the limiting size distribution CAr(n) of the leader for large N. Integrating Eq. (IHT 
over V and writing z = e^'^ we get 



nje 
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where 
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dv'P{n,v) exp[— y(s,z;)] 



duV{k, u) 



(32) 



(33) 



with 'P{k, u) defined in Eq. Now, for large N, we need to investigate the behavior of the r.h.s of Eq. (|32|) in the 
limit s ^ 0. It is evident that to obtain a sensible limit of the r.h.s in Eq. ([5^ one needs to consider the scaling 
limit s ~f but n oo keeping the product sn finite. In this scaling limit, it is not difficult to see that the dominant 
contribution to the sum Y(s, v) in Eq. p3p comes from those terms where the product sk is finite in the limit s — > 0, 
i.e., terms where fc 1/s is large. 

Now, for large k, it follows from the definition in Eq. ^ that for a symmetric (/>(w) with zero mean and a finite 
variance cr^ = v'^(f>{v)dv is finite, one can invoke the central limit theorem to assert that V{k, u) has a Gaussian 
distribution, i.e.. 



V{k,u) 



exp 



kv? 
2^ 



(34) 



Substituting this result in Eq. (133p and replacing the resulting sum by an integral in the small s limit, we find that 
an appropriate scaling limit of Y{s, v) exists when s — > 0, w — > keeping the ratio v/yG fixed. In this scaling limit, 
we get 



dy 



rfc 



The integral on the r.h.s of Eq. (|35p can be done explicitly and one gets to leading order in small s 



Y{s, v) « -2*-^ ^ ^ ^ + + 1 



(35) 



(36) 



where C = 0.57721566 ... is the Euler's constant and w — is the scaling variable. Finally we substitute this 

expression for Y{s, v) and the limiting Gaussian form of V{n, v) from Eq. onto the r.h.s of Eq. ((5^ . perform the 
subsequent integration and finally arrive at the following expression in the scaling limit s — > 0, n — > oo but keeping 
the product ns fixed 



VC„(n)e-^^«^G(ns) 



where the scaling function G(.x) is given by 

G(x) 



2b 6-=^ 
x/tF x/x 



dye" 



X 



(37) 



(38) 



and 



& = 6*^/2 = 1.33456... . (39) 

We notice one remarkable fact: even though we needed to have finite in arriving at the result in Eq. cr^ has 

dropped out of the final expression in Eq. (j38p . Thus for any finite ct^, the scaling function G{x) is universal and is 
actually independent of the value of . 

From the expression in Eq. psp it is easy to compute the moments of the leader size distribution for large N . For 
example, to compute the fc-th moment {n^') for large iV, we multiply both sides of Eq. (I38p by and sum over all 
n. This gives, 

^(n^Ve-^^«^E^"'""'"W dye-y' ^l + y^/ns. (40) 
V n Jo 

We next make a change of variable y = y^sn u in the integral and carry out the sum over n. In the scaling limit, this 
sum can be replaced by an integral which can be explicitly performed and we get for small s 

^ , .AT,. 26r(fc + i/2) r du _ r(fc-i/2)r(fc + i/2) 

^^^^ ^ V^s'=+V2 {l + u^)k-^ Yik) ' • ^^'^ 

Note that this expression is strictly valid for fc > 1. For k = 0, one can show independently that J2n C'Ar(n)e~*^ ^ 1/s 
as s ^ as expected since the PDF CN{n) is normalized to unity. Inverting the Laplace transform in Eq. (|¥T|l .it 
then follows that for large TV, the fc-th moment of the leader size behaves for fc > 1, 

A, 7^=^1/2. ^here A, = ^Iti/^. (42) 

r(fc) 

In particular the average leader size (fc = 1) behaves for large N as 

{n)^ « b^/^N^/^ = (2.63533...)iVi/2 (43) 

From the expression for the moments in Eq. (|4ip . it follows that for large iV, the leader size distribution C]\[(ji) 
has the following scaling form 

(44) 

where the scaling function W{x), for < a; < 1 is such that, for all fc > 1, 

W{x)x''dx = b^^^^^—\^ (45) 

r(fc) 

Finally, one can analytically continue the Eq. (|45p for noninteger fc and then invert the equation to obtain the 
following explicit expression of W{x) valid for <C .t < 1 

W{x)^-^x-^^^{l-x)-^/^ (46) 

It is easy to verify that for all fc > 1, the above expression for W{x) does indeed satisfy Eq. (|45p . Note that the scaling 
in Eq. (|44p is valid only for 1 <^ n N. In particular it does not hold when n ~ 0(1)- This is also evident from the 
explicit form of the scaling function W{x) in Eq. (|46p which diverges as x~^^'^ as x — > 0, making it non-normalizable. 
Indeed, for small n, one has to keep the full expression of the generating function to recover the correct normalization. 

Thus the main result of this subsection is the universal scaling form of the leader size distribution in Eqs. (|44p. 
()45p and (I46p valid for any continuous and symmetric (j){v) with a finite cr^ and remarkably it is independent of the 
actual value of the cr^ . 

B. The Limiting Leader Velocity Distribution 

The marginal velocity PDF of the leader Diy{v) can be obtained by summing Eq. (|3ip over n. This gives 



°° ^ 1 If 

Dn{v) = -— p^{v) exp - / 

N=l ^ ^ L 



Pz{u)du 



(47) 



A somewhat more convenient quantity is the cumulative velocity distribution of the leader, Fn{v) 
One can easily derive its generating function by integrating Eq. (|47|) 



N=l 



Fn{v)7J 



1 - z 



f_^D^{v')dv'. 



(48) 



Note that as u — > cx), the r.h.s of Eq. (|48p becomes z/{l — z). This is consistent with the l.h.s of Eq. ((48|) since 
Fn{v) — *■ 1 as w — *■ oo and thus the generating function in that limit is precisely z/{l~ z). 

The result for the cumulative velocity distribution of the leader in Eq. (|48p is exact for all N . We next address the 
question of the limiting leader velocity distribution as N ~* oo. For large N , one needs to investigate the r.h.s of Eq. 
()48|1 in the limit z — > 1. In that limit, the r.h.s scales as exp [— pi{u)d'u\ /(I — z), provided the numerator is finite. 
In such cases, it follows from Eq. that Fn{v) tends to a limiting A^-independent distribution Foo{v) given by 



Foo{v) = exp 



pi{u) du 



exp 



oo -. 

^ n 

n=l 



Vin, u) du 



(49) 



Note that for some (f>{v) this limiting distribution may not exist. In fact, we will see later that when (f>{v) has the 
Cauchy distribution, which is symmetric and continuous, there is no iV-independent limiting distribution. It is evident 
from Eq. (jH]) that, unlike the limiting leader size distribution, the limiting leader velocity distribution, whenever it 
exists, is highly nonuniversal and depends explicitly on the initial velocity distribution (/)(w). 

To see how this limiting distribution in Eq. (j49p looks like, we first consider the Gaussian distribution, 0(u) = 

exp(— w^/2CT^)/\/27rCT^. In this case, V{n,u) — y/^^j exp — exactly for all u. Substituting on the r.h.s of Eq. 

()49p and performing the integral one can express the limiting distribution Fadv) as a function of the dimensionless 
variable v/a, 



F^iv) ^ Fiv/a); 



where F{z) — exp 



(50) 



While we were unable to perform the sum exactly in Eq. ([50|) . the function F{z) can be easily plotted using Mathe- 
matica and is shown in Fig. It is easy to show from Eq. (|5n|) that the function F{z) has finite nonzero support 
only for positive z. It is identically zero for all z < 0. Furthermore, the asymptotics of F{z) as z — > and z ^ oo 
can also be derived. We find 



F{z) w \/2i 



2ttz 



as z 

as z ■ 





oo. 



(51) 
(52) 



The asymptotic behavior as z ^ oo is easy to derive, as in this limit the dominant contribution comes from the n = 1 
term of the sum in Eq. ([50)1 . In contrast, the other limit z ^ in Ec^. (|5ip is more tricky to derive. We provide a 
derivation in Appendix-B. 

Let us make a couple of interesting observations on the general formula for the limiting velocity distribution in Eq. 
(l49ll when it exists. 



(i) For any continuous and symmetric initial velocity distribution 0(w) for which the limiting velocity distribution 
Foo{v) exists, Foo{v) = for any v < 0. This follows from the fact for w = 0, by symmetry, V{n,u) du — 1/2. 
Thus, the sum in Eq. (|49p diverges for v ~ and -Foo(O) = 0. Since Fac{v) is a cumulative distribution, it is a 
nondecreasing function of v. Hence it follows that if Foc{v = 0) = then Foo{v) = for all v < 0. The implication of 
this result is rather interesting. It implies that even though initial velocity distribution may be symmetric in v and 
there may be a lot of particles initially with a negative velocity, the eventual velocity of the leader, i.e., the rightmost 
cluster is always positive. 

(ii) Note that the velocity of the leader is also the maximum of the final velocities. Now if the collisions were elastic, 
the particles would have merely interchanged the velocities in each collision (this is called the Jepsen gas), and in the 
final state the velocity of the leader would have been the maximum of all the initial velocities [2^. Therefore, in this 
case the distribution of the leader velocity would have been given by the usual extreme value statistics of uncorrelated 
random variables The velocity of the final leader of elastic Jepsen gas scales as u « w + for large A^, 
and the limiting distribution of the scaled velocity w has only three possible forms depending on the tail of (j){v) [26l |. 
However, the scaling parameters a^r and depends on N as well as on the complete form of 4>{v). For example. 
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FIG. 4: (Color online). The limiting cumulative velocity distribution Foo{v) = F{z = v/a) plotted against z, for the Gaussian 
distribution ^(u) = exp{—v'^/2a^)/\/2Tva^. 



when (j){v) is Gaussian, gn = {2lnNy^/^ and In = (21n7V)i/2_ (21n7V)"i/2(-inij^^_l_lj.^4^)/2. The limiting velocity 
distribution is Gumbel in this case. Thus, unlike the elastic gas where the leader velocity increases with N, it follows 
from Eq. I|49p (when the limiting distribution exists) that in the sticky gas, the leader velocity remains of 0(1) in the 
large N limit. This is due to the strong correlations generated between the velocities of different clusters in the fan 
state. Thus Eq. (|49|) provides an exactly solvable case of extreme value statistics of correlated random variables. 



Exact Leader Statistics for the Cauchy Distribution 



In this subsection, we consider the special case of the symmetric Cauchy distribution for the initial velocity. 



1 



a 



(53) 



For this distribution, the mean is zero, but the second moment diverges. So, the results of the previous two 
subsections, where it was assumed is finite, do not hold. However, the principal result in Eq. pip is still valid and 
we show here that the Cauchy distribution preents an excatly solvable case in the sense that the r.h.s of Eq. (PT|) can 
be exactly evaluated in closed form. 

The first simplification occurs when one calculates V{n,u) defined in Eq. (jj]). Since the Cauchy distribution is 
stable, it turns out (as can be easily proved) that 



1 



V{n, u) — (j){u) — — 



TT + 



(54) 



Thus P{n, u) is completely independent of n. This particular feature of the Cauchy distribution was also used recently 
to obtain exeat results for the statistics of records in a sequence of random walk in presence of a drift [1^ . Using this 
independence on n, it then follows from the definition in Eq. ([7|) that 



Pziu) 



One can then trivially do the integral 



TT- — o Infl — z). 



arctan(u/a) ln(l — z). 

2 TT 



(55) 



(56) 



Substituting these results in Eq. (|3ip gives the foUowing exact expression 



1 a z" 1 1 

LN(n, v)z^ ^ — (1 - where ^liv) = - + - arctan(w/a). (57) 

^ — ' TT + n 2 TT 

N=l 

Comparing coefRcient of gives us the exact joint distribution of size and the velocity of the leader 

LN[n, v) = - — — — , ..-p. , yy (58) 

Marginal size distribution: To compute the marginal size distribution CN{n) = Lpf{n,v) dv of the leader, it 
is convenient to do the integration over v in Eq. (|57p . This gives the exact generating function 

°° n+l 

J2CMin)z''^ \ ^. (59) 

^-^ nil — z) m 1 — z) 

From this generating function, it is easy to calculate the moments {n^) of the leader size for large N by investing 
the singularity of the r.h.s in Eq. ([59|) near z = 1. Skipping the details, we find that for large N and for fc > 1 

This result should be compared to the case in Eq. where is finite. In particular, the average leader size grows 
as {n) « N/\nN, much faster than the N^^^ growth for the case where is finite. 

Marginal velocity distribution: Similarly, the marginal velocity PDF of the leader Dm{v) = X]„ L]^{n, v) can be 
computed conveniently directly from Eq. ((57|) . We get 



^ D^{v) z^^- [- ln(l - z)] (1 - z)-''W, (61) 

Af=l 

from which one can compute the generating function for the cumulative velocity distribution Fn{v) = Dn{v') dv' . 
We get a nice compact expression 

C30 



^F7v(i;)z^ = (l-z)-^(")-l. (62) 

N=l 

Comparing powers of z^ gives a very simple but nontrivial distribution valid for all A'^ > 1 

, , T(u.(v) +N) , ,,11 , , , 

Fn{v) = Y^j^^^-j ; where fi{v) = ^ + - arctan(w/a). (63) 

Note that when v oo, fj,{v) —f 1 and when v — oo, iJ,{v) 0. Thus, in the limit v —^ oo, F^{v) 1 and 
when V — oo, Fn{v) ^ as expected. Also, for = 1, we get Fi{v) = ix{v) which is just the cumulative Cauchy 
distribution. This is expected, because the velocity of a single particle remains unchanged as there is no collision. 

We note from this exact velocity distribution in Eq. ((55)) that the distribution is TV-dependent explicitly and does 
not have any nontrivial iV-independent limiting form as iV — > oo, unlike in the previous subsection for the Gaussian 
case. In fact, by ananlysing the singularity near z = 1 in Eq. ()62p . we see that for any fixed as N oo, Fi^{v) 
decays with A'^ as a power law with an exponent that depends continuously on v 

Since the exponent 1 — ^{v) ~ 1 / 2 ~ avctaniv / a) / n decreases monotonically from 1 (as v — s- — oo) to zero (as v — > oo), 
it follows that the distribution for smaller values of v tends to zero faster than the distribution for higher v. Thus the 
marginal leader velocity distribution in the Cauchy case is very different from that of its Gaussian counterpart. 



V. DISTRIBUTION OF THE TOTAL ENERGY IN THE FAN STATE 



In this section we compute the distribution of the total energy of the clusters in the fan state. For simplicity, we 
restrict ourselves to the initial Gaussian distribution of the velocities, (f>{v) = exp(— i;^/2)/\/27r where we have also 
set cr^ = 1 without any loss of generality. Thus, initially, the total (kinetic) energy of the system, Eq 
distributed according to the PDF 



Pn{Eo) 



E. 



N/2- 



T{N/2) 



(65) 



Now, as the system evolves, it dissipates energy through collisions, and once the fan state is reached there is no more 
energy loss. It is natural to ask how the PDF of the energy evolves with time from its initial form in Eq. (j65p . and 
in particular how does this energy PDF look like in the final fan state beyond which there is no more dissipation. 
Let Qm{E) denote the PDF of the total energy in the fan state, i.e.. 



Qn{E) 



5\E 



1 ^ 



(66) 



where the angle brackets (■ • • ) denote the average over the random variables k, {ni} and {ui}, whose joint PDF is 
given by Eq. It is convenient to consider the Laplace transform of the energy in Eq. ([55]) 
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exp 
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i=l 



(67) 



To evaulate this average using the measure in Eq. ([2]), we need to first evaluate 'P{n,u) in Eq. For Gaussian 



4>{v), this is simple: V{n,u) 



exp 



n u 
2 



exactly for all n and u. We substitute this in the measure in 



Eq. Q and note that calculating the average in Eq. (j67p just amounts to renormalizing the effective V{n,u) 
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n u 
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With this input, one can then carry out exactly the same procedure as in subsection III A 



and one arrives at the final result 
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(68) 



where we recall Si{N, k)/Nl from Eq. (114p . The inverse Laplace transform yields 
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Qn{E)^Y. 







'Si{N,k)' 
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m 



(69) 



Since in the final state the number of clusters k is distributed according to pN(k) — Si{N, k)/Nl, it is interesting to 
note that, for a given number k of final clusters, the distribution of the total energy is identical to that of k initial 
particles. 

Equation (|69p provides an exact expression for Qiy{E), for all values of E and N. However, it is always desirable 
to have an closed-form expression, even though its validity would require taking the limit of large N. To achieve this 
goal, we first express Eq. in the following form by using the generating function given in Eq. 



)^iE)e--^dE = l^^^ 



r i/vr 



a 



(70) 



At this point, we make a short detour to compute the mean and variance of the final total energy, as they are 
readily available from the first and second derivatives of Eq. ((70|) with respect to a at a = 0. We find the mean as 



(E) 



N 



(71) 



where Hiq = X]fc=i ^ ^ is the harmonic number, and asymptotically, Hj^ ^ IniV + C where C is the Euler's constant. 
The variance is 

{E^) {Ef = (72) 

where ^ — X^a^i ^ ' generalized harmonic number. Note that Hn = Hn,i- The limiting value of H^.r as 
oo is the Riemann zeta function i/oo,r — C{^)i and C(2) = 7r^/6. 
We now return to Eq. (ffO)) . and proceed along the main course to obtain a closed-form expression of Qn{E) for 
large N . The inverse Laplace transform of Eq. (|70p is formally given by the Bromwich integral 

Qn{E)^— e^("Ma, (73) 

where 

5(q) =lnr(^iV + 1/VTT^) -ln(Ar!) 



- In r (^1 / VTT^ j + aE. (74) 
Now a saddle point approximation of the integral in Eq. (j73p yields 



Qw(i?)«e^("*)/v/27r|^"(a*)| , (75) 

in which a* is determined implicitly in terms of E and by the saddle point condition S'{a*) = 0. Enforcing this 
condition to Eq. ([7i)) gives the saddle point equation for a* as 



2 



^(7V + /3) - V(/3) , (76) 



where /? = l/Vl + a* and tp{x) — [d/ Ax)\nT{x) — r'(a;)/r(x) is the digamma Function. As %1){N + /3) ~ InA^, 
Eq. ([7H|) suggests that when both E and are large, the natural scaling variable is z = 2£'/lnA^. In terms of this 
scaling variable, for large N ^ from Eq. (|76p we find 

Now, employing Stirling's approximation (for large N) in Eq. (|74p and using Eq. (j77p . in terms of z we have 



5(a* 

Similarly, we find 



3 1 z 

77^ ^ -77-1 
2 2 



lniV-lnr(zV3) + o(J^) 



(78) 



^"(a*) - ^z^/MnTV. (79) 



Therefore, substituting S and S" form Eq. ((75|) and Eq. ([7^ respectively, in Eq. ([75]) . finally we obtain 



-lnA.j.(^-lnivj ^,7,^, (80) 



where the large deviation function h{z) is given by 



h{z)^l+'--\z''\ (81) 

The saddle-point approximation Eq. (|80p is indeed a very good for large N , as we demonstrate in Fig. [5] by comparing 
with the exact result Eq. ^ for N = 2000. 

Near z = 1 we have h{l + x) — a;^/6 -|- 0{x^). Therefore, from Eq. ((501) '^^ find that Qn{E) has a Gaussian form 
close to the mean {E) ^ (1/2) In TV, with a variance (i?^) — (E)-^ ~ (3/4) In A^; — we note that by taking the large A^ 
limit in Eq. (I7ip and Eq. (|72p respectively, we independently recover these mean and variance. However, far away 
from the mean, the PDF is quite asymmetric as clearly seen in Fig. \5\ 
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FIG. 5: (Color online). Main: Distribution of the total energy for A'^ = 2000. The solid (red) line plots the saddle-point 
approximation Eq. (|80|l . and the dashed (blue) line plots the exact form Eq. (|69p . evaluated using Mathematica. Inset: 
The solid (red) line plots the large deviation function Eq. (|8ip . and the dashed (blue) line plots — (In 7V)~^ In Qjv (| In Af) — 
(IniV)"^ ln[(T" ^'O-Nf^'^z^''^ ^{z^'^)] as a function of z, using the exact form of Qn(E) given in Eq. 

VI. SUMMARY AND CONCLUSION 

In summary, we have studied analytically various statistical properties in the final fan state of a one dimensional 
sticky gas of N particles. The particles are initially distributed randomly in space, each with unit mass and with an 
initial velocity drawn independently from an identical distribution (^(w). Each particle moves ballistically and when 
two particles collide, they form a single cluster with the total mass and total momentum conserved in the collision 
process. In the long time limit, the system reaches a fan state where the system consists of a finite number of clusters 
with their velocities increasing from left to right and there is no further collision. We have shown that the sizes and the 
number of clusters are distributed universally in the fan state, independent of <f>{v) and this distribution is identical 
to that of the cycles lengths in random permutation problem. We have also computed exactly the distribution of the 
size and the velocity of the rightmost cluster (leader) moving with the largest velocity. Furthermore, the distribution 
of the total energy in the fan state is also computed exactly for Gaussian </'(w). Our results provide an exactly solvable 
case of a many-body system (of finite size N) with dissipative dynamics and in particular, brings out the universal 
features of the final stationary state in an explicit way. 

There are interesting open issues for further research. Here we have focused solely on the final fan state. It would 
be interesting to study the dynamics, i.e., the approach to this fan state and investigate to what extent the universal 
features exist away from the fan state. 

There are also immediate generalizations of the problem studied here that are open for future research. For example, 
it would be interesting to know to what extent the universal features in the fan state remain valid for an inhomogeneous 
sticky gas. The inhomogeneity can arise either from unequal initial masses or nonidentical velocity distributions in 
the initial condition. For an infinite system, i.e., with a finite density of initial particles, the ballistic aggregation 
model with inhomogeneous initial condition has been studied [11] . It would be interesting to extend this study to the 
case of a finite number N of particles. 

Finally, it would be interesting to study this finite system of N particles when the collisions are inelastic, but not 
necessarily sticky, i.e., with a nonzero coefficient of restitution. This problem has been studied for an infinite system 
and many asymptotic properties such as the energy decay and the velocity distribution were found to be similar to 
the sticky gas limit [29, 30]. It would be interesting to see to what extent the universal features in the stationary 
state for finite N are retained when one introduces a nonzero coefficient of restitution. 
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APPENDIX A: PROOF OF EQUATION (fTsj) 



In this Appendix, we derive the distribution of the cluster sizes (regardless of their positions) starting from joint 
probability distribution of the cluster sizes and velocities given in Eq. To calculate Pr{ci, C2, . . . , cjv} from the 
joint-PDF pN^k] {ni,Ui}), we must: (i) integrate out the velocities, (ii) sum over all configurations with the same 
cluster sizes (but placed in different orders) (iii) change the variables from the list (ni, , . . . ,nk) that represents the 
sizes of the k consecutive clusters to the set {ci, . . . , cn} that encodes the number of clusters of a given size. We thus 
have 

Pr{ci,C2, . . . ,C7v} = ^ pN{ni,...,nk), (Al) 

Permutations of 
unequal size clusters 



where the sum is restricted to permutations that exchange clusters having different sizes i.e. permutations of the list 



(ni, . . . ,nk) that exchange n^'s having different values. Using Eq. ^ and the fact that Y\i=i = Il^i J'^^ i find 
that pAr(ni, . . . , rifc) is given by 

PN{ni, ■ • ■ ,^^fe) 

From the definition of V{ni^ui) given in Eq. (j4]), the integral on the r.h.s of Eq. (|A2p is found to be 

„ Af k-l I Wf+i Nt \ 

\{dv.^{v.)\{e[—- E -A (A3) 
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^ ^ / n duiVim, m) Jl e{ui+i - ue) . (A2) 
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where Ni — ni + . . . + ni for £ — 1, . . . , /c. The proof of Eq. psp hence reduces to showing the following identity: 

N k-l I A^f+i Nit 
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unequal clusters 
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ci!c2!...CAr! ' 



(A4) 



where we recall that c\ is equal to the total number of I's in the list (ni, . . . , nfc), ci is equal to the total number of 
2's in this list etc... This identity is proved by the following argument: consider a permutation that exchanges two 
clusters of the same size p; then, by relabeling the Vi (which are dummy variables) we observe that the integral on the 
left hand side (l.h.s) of Eq. (|A4|) is invariant. There are Cp! possible permutations of the Cp clusters of size p. Thus 
we have formally 



cx\ci\ . . .Cat! 



Permutations of All Permutations 
unequal clusters 

If we have k distinct numbers mi . . . Ufe there is only one permutation that orders them in increasing order. Therefore, 
we have 

k~l 

E n ^("'^(^+1) - '"'t(£)) = 1 , (A6) 

o-eEfc 1=1 

where is the symmetric group of k elements. After substituting this identity in Eq. (|A5|) and using the fact that 
the PDF (\>{yi) is normalized, Eq. (IA4|) is finally proved. We note that the method of proof used here is closely related 
to the derivation by F. Spitzer of a generalization of the Sparre- Andersen theorem [15| . 



APPENDIX B: LIMITING BEHAVIOR OF F{z) AS 2 IN EQ. ^ 
Let us write F{z) in Eq. (|50p as F{z) = exp[— ^(z)] where the sum 



n—1 ^ 



We want to compute the behavior of S{z) as z ^ 0. We first note the identity 

1 - e-^'/2" 



Subtracting Eq. from Eq. (|BT|) gives 



Siz) + - In 



1 - e 



y- 

^ 2n 



erfc I A / 2 ^ 



(Bl) 



(B2) 



(B3) 



Now, let us define a new variable y = nz'^/2. As n changes by 1, y changes by Ay = z'^/2. Now, in the limit z ^ 0, 
this increment Ay is very small. Thus, one can replace the sum over n on the r.h.s of Eq. (jB3|) by an integral over y 
to leading order for small z, giving 









rdy 








Jo y 



[erfc(V^)-e-^] 



The integral / on the r.h.s is just a constant and can be evaluated by parts 



1 



/ e-yHy)dy-^ 
Jo Vi" Jo 



My) 



e ^ dy 



(B4) 



(B5) 



The two integrals on the r.h.s are elementary ones and can be easily evaluated to give finally, I — ~ ln(2). Substituting 
this in Eq. jBSj and taking the z — > limit gives 



implying as z ^ to leading order 



S{z) -> -ln(\/2z), 
F{z) = exp[-S'(z)] -> V2; 



(B6) 



(B7) 
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